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I. INTRODUCTION 



In the past few years it has become clear that topological defects produced in the early 
universe may have a considerably richer microstructure than had previously been imag- 
ined [Q. In particular, the core of a defect acquires additional features at each subsequent 
symmetry breaking which preserves the topology of the object. The new microphysics as- 
sociated with additional core structure has been exploited by several authors to provide a 
new, defect-based scenario for electroweak baryogenesis 0,0. 

The purpose of the present paper is to constrain general particle physics theories by 
demanding that the microphysics of defects in these models be consistent with the require- 
ments of the standard cosmology. The basic idea, due originally to Davis and Shellard , 
is as follows. If a spontaneously broken field theory admits linear topological defects - cos- 
mic strings - which subsequently become superconducting, then an initially weak current 
on a closed string loop will automatically tend to amplify as the loop undergoes dissipative 
contraction. This current may become sufficiently strong to modify the dynamics and halt 
the contraction so that the loop settles down in an equilibrium state known as a vorton. 

The population of vorton states produced by such a mechanism is tightly constrained by 
empirical cosmological considerations. It was first pointed out by Davis and Shellard that 
to avoid obtaining a present day cosmological closure factor i7 greatly exceeding unity, any 
theory giving rise to stable vorton creation by superconductivity that sets in during string 
formation is ruled out if the symmetry breaking scale is above some critical value. One of the 
first attempts to estimate this critical scale indicated that it probably could not exceed 
that of electroweak symmetry breaking at about 10^ GeV by more than a few orders of 
magnitude. Such strong limits are of course dependent on the supposition that the vortons 
are absolutely stable on timescales as long as the present age of the universe. However, 
even if the vortons only survive for a few minutes, this would be sufficient to significantly 
affect primordial nucleosynthesis and hence provide limits of a weaker but nevertheless still 
interesting kind 0. 
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In all previous work it was supposed that the relevant superconductivity sets in during 
or very soon after the primary phase transition in which the strings are formed. What is 
new in the present work is an examination of the extent to which the limits discussed above 
are weakened if it is supposed that superconductivity sets in during a distinct secondary 
phase transition occurring at what may be a very much lower temperature than the string 
formation scale 

The structure of the paper is as follows. In Section IIA we shall, for completeness, 
give a brief introductory review of string superconductivity. In Section IIB we describe 
the mechanism of formation of a vorton from an originally distended string loop and in 
Section IIC we summarize the basic properties of vorton equilibrium states. In Section 
IIIA we first comment on how it can be that the formation scale and the superconductivity 
scale can be separated by many orders of magnitude. We then demonstrate, using a suitably 
simplified statistical description of the string network, how to estimate the vorton abundance 
for a generic theory as a function of the temperature and the symmetry breaking scales. In 
Section IIIB we apply this procedure to the relatively simple case when the superconductivity 
develops during the early period when dissipation is mainly due to the friction of the ambient 
medium. In Section IIIC we go on to treat the more complicated situation that arises if the 
superconductivity does not develop until the much later stage in which dissipation is mainly 
due to gravitational radiation and in Section HID we briefly comment on stability issues. In 
Section IV we consider the comparitively weak bounds that are obtained if it is supposed 
that the vortons are stable only for a few minutes. Finally, in Section V we consider the 
rather stronger bounds that are obtained if the vortons are of a kind that is sufficiently 
stable to survive as a constituent of the dark matter in the universe at the present epoch. 
We conclude in Section VI. 

II. CONSEQUENCES OF COSMIC STRING SUPERCONDUCTIVITY 
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A. Currents in the Witten model. 



In so far as it has been developed at the present time, the quantitative theory of vorton 
structure has been entirely based on the supposition that the essential features are describ- 
able in terms of a simple bosonic superconductivity model of the kind introduced by Witten 
|T^. This category of models consists of spontaneously broken gauged U{1) xU (1) field the- 
ories, which generalise the even simpler category of spontaneously broken gauged U{1) field 
theories on which the standard Kibble description of non superconducting cosmic strings is 
based. 

The Kibble model is characterised by a potential V with a quartic dependence on a 
complex Higgs field of the familiar form A(|</)p — v'^y. Here A is a dimensionless coupling 
and f is a mass scale of order the "Kibble mass" m^, whose square is identifiable with the 
string tension, T, which is (in this model) constant and equal to the mass per unit length. 
The string is defined as the region in which |0| is topologically excluded from its vacuum 
value as given by ^ m^. 

In addition to 0, the Witten model contains a second complex scalar field, a, and is 
characterised by a quartic potential depending on several dimensionless parameters that, 
like A, are assumed to be of order unity. Witten's potential function also depends on a 
second mass parameter, m„ say, which determines the temperature scale below which a 
gives rise to a current carrying condensate on the vortex. In the Witten model the vortex 
defects are cosmic strings in which the tension T is no longer constant but variable, as a 
function of the current magnitude \j \ , attaining its maximum (Kibble) value only when the 
current magnitude vanishes, so that in general one has T < m^. 

In earlier discussions of vorton physics it was implicitly or explicitly supposed that 
the magnitude of rricr was not very different from that of the original symmetry breaking mass 
parameter m^. In that case, the formation of the cr condensate could be considered as part of 
the same symmetry breaking phase transition as that by which the strings themselves were 
formed. Our purpose here is to consider scenarios in which rricr may be very much smaller 
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than mx, as will occur when successive phase transitions at two entirely distinct cosmological 
epochs are involved 0. It is only after the second phase transition that a condensate with 
amplitude \a\ and angular variable phase, 6 say, will form on the string world sheet. There 
then exists an identically conserved worldsheet phase current with components 

r = ^^"'dbe , (1) 

where e"''^ are the components of the antisymmetric unit surface element tensor induced on 
the 2-dimensional worldsheet. As well as this topologically conserved phase current, there 
will also be a dynamically conserved particle number current with components given |0 by 

ja = 2t{dae - eAa) , (2) 

where S is the surface integral of |o"p over the vortex core cross section. Here Aa are 
the induced components of the electomagnetic background potential and e is the coupling 
constant associated with the carrier field if it is gauged. When e is non zero will be gauge 
dependent, but is physically well defined and will determine a corresponding electric 
surface current given by 

r = ef . (3) 



B. Formation of Vorton States. 

In the case of a closed string loop the conserved surface currents characterised above will 
determine a corresponding pair of integral quantum numbers that are expressible in terms 
of circuit integration round the loop. These are given by 

N = frdL, Z = jfMa. (4) 

where Ma are the components of the length element normal to the circuit in the worldsheet. 
Note that even when j'^ is gauge dependent, N is well defined. A non conducting Kibble 
type string loop must ultimately decay by radiative and frictional drag processes until it 
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disappears completely. However, since a Witten type conducting string loop is characterised 
by the classically conserved quantum numbers and Z, such a loop may be saved from 
disappearance by reaching a state in which the energy attains a minimum for given non zero 
values of these numbers. 

In view of a widespread misunderstanding about this point, it is to be emphasised that 
the existence of such energy minimising "vorton" states does not require that the carrier 
field be gauge coupled. If there is indeed a non vanishing charge coupling then the loop will 
of course be characterised by a corresponding total electric charge 

Q=<fPdia, (5) 

in terms of which the particle number will be expressible directly as Z = Q / e. However, the 
important point is that even in the uncoupled case, for which and hence also Q, vanish, 
the quantum number Z will nevertheless remain perfectly well defined. 

Although the essential physical properties of a vorton state will be fully determined by 
the specification of the relevant pair of integers and Z, it is not true that any arbitrary 
choice of these two numbers will characterise a viable vorton configuration. This is because 
the requirement that the strictly classical string description should remain valid turns out 
to be rather restrictive. To start with, it is evident that to avoid decaying completely like a 
non conducting loop, a conducting loop must have a non zero value for at least one of the 
numbers and Z. In fact, one would expect that both these numbers should be reasonably 
large compared with unity to diminish the likelihood of quantum decay by barrier tunneling. 
However, even for moderately large values of and Z there will be further restrictions on 
the admissible values of their ratio Z jN due to the necessity of avoiding spontaneous particle 
emission as a result of current saturation. 

The existence of a maximum amplitude for the string current was originally predicted 
by Witten himself []T0[ . However, quantitative knowledge about the current saturation phe- 



nomenon remained undeveloped until the appearance of an important pioneering investiga- 
tion by Babul, Piran, and Spergel |]l2l who undertook a detailed numerical analysis of the 



mechanism whereby the presence of a current tends to diminish the string tension T. In a 
non conducting string, the tension T is identifiable with the energy per unit length £^ but 
in the conducting case the diminution of T is accompanied by an augmentation of £ such 
that 

T <ml<£ . (6) 



The analysis of |T2| provided an empirical "equation of state" specifying T as a non-linear 
function of the current magnitude \j\ and hence of the energy density £. The minimum of 
T, and hence the maximum allowed value for £, is obtained when the current amplitude |j| 
reaches a critical saturation value with order of magnitude 

\3\'^£-r<ml. (7) 

Such knowledge of the equation of state is precisely what is required for investigating string 
dynamics. 

We now apply this to vorton equilibrium states. The phase angle, 9, is expressible in 
terms of the background time coordinate t and a coordinate ^ representing arc length round 
the loop by 

e = ut-ki , (8) 

where lo and k are constants. If the total circumference of the vorton configuration is denoted 
by £v then it can be seen that these constants will be given in terms of the corresponding 
quantum numbers by 

uj = ^— , k = — — . (9) 
The general condition [jl4| for purely centrifugally supported equilibrium is that the rotation 



velocity should be given by a formula of the same simple form as the one |jT3[ for the speed 
of extrinsic wiggle propagation, namely 

v' = j. (10) 

Note that this relationship applies when the electrical coupling is absent or, as will commonly 
be the case, negligible because of the smallness of the fine structure constant [|l^ . 
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C. Properties of Vortons. 



The random distribution of initial values of the quantum numbers N and Z leads to 
the formation of a range of qualitatively different kinds of vorton states. We shall now 
briefly review the reasons why it is to be expected that the most numerous initially will be 
of approximately chiral type, meaning that their rotation velocity is comparable with the 
speed of light as given hj v = 1. Note however, that it may be that vortons of the very 
different subsonic type, with v <^ 1, will be favoured by natural selection in the long run. 

In the case of a spacelike current, with u"^ < k"^ (the only possibility envisaged in [0), 
the velocity given by (|I0[) is to be interpreted simply as the phase speed: 

uj -kZ 

V = — = ^ — . (11) 

k EN ^ ^ 

Since we are assuming that rricr is small compared with rrix, the saturation limit (|^) implies 

^ « 1 , (12) 

which obviously, by ([TIID, implies the property of approximate chirality, w ~ 1. It is expected 
on dimensional grounds that, although the sectional integral E is a function of the current, 
it will not get extremely far from the order unity. This behaviour has been observed nu- 
merically in particular cases (see also |T^JT7[] ). It can therefore be deduced that 



approximate chirality requires the vorton to be characterised by a pair of quantum numbers 
having roughly comparable orders of magnitude, 

\Z\^N, (13) 

where we use an orientation convention such that is always positive. 

On purely statistical grounds one would expect that the two quantum numbers would 
most commonly be formed with comparable order of magnitude, so as to satisfy (|T3]), though 
not necessarily with an exact ratio small enough to give a spacelike vorton current. In fact, 
although the limit definitely excludes the possibility of vorton states with ^ \Z\, 
it turns out that there is nothing to prevent the existence of non-chiral vortons having a 



current that is not just marginally timelike (c<j^ > fc^), but which have \Z\ ^ A^. In the 
timelike case the velocity v in (p!0D is not the phase speed but the current velocity given by 

v = k/uj . (14) 

A preliminary application |^0| of the principles described in the following sections sug- 



gests that subsonic vortons, which will necessarily be characterised by |Z| ^ A^, will initially 
be formed in much smaller numbers than the more familiar chiral variety described by (|13[). 
This means that if ordinary chiral vortons are sufficiently stable to survive over cosmologi- 
cally significant timescales (a question that must be left open for future research) then they 
can provide us with much stronger constraints on admissible particle theories than the more 
exotic subsonic variety. In order to keep our discussion as clear and simple as possible, we 
shall therefore say no more about subsonic vortons and restrict our attention to the chiral 
variety. 



Chiral vortons are only marginally affected by the electromagnetic coupling, e [[15 
Therefore, they can be described by the simplest kind of elastic string formalism in which 
electromagnetic effects are ignored altogether. This means [|1^] that the mass energy of 
such a vorton will be given in terms of its circumference by 

= 4(^ + T) . (15) 

In order to evaluate this quantity all that remains is to work out l^. In the earliest work 
on vorton states it was always presumed that they would be circular, with radius therefore 
given by = £v/27r and angular momentum quantum number J given [|l^] by J = NZ or 
equivalently by = STi^/An'^. Thus, eliminating J, one obtains the required result as 

4 = (27r)i/2|^Z|i/2^^^)-i/4 ^ (27r)i/2|Arz|i/2^-i . (16) 



More recent work |^T| has established that even in cases where the vorton configuration is 
strongly distorted, the expression ( [TB| ) will remain perfectly valid. Combining this with (|I3|), 
and recalling that for chiral vortons \Z\ ^ N, we thus obtain a final estimate of the vorton 



mass energy as 

8 



The preceding formulae are based on a classical description of the string dynamics. This 
is valid only if the length is large compared with the relevant quantum wavelengths, of 
which the longest is the Compton wavelength associated with the carrier mass mo-. It can 
be seen from (0) that this condition, namely 

4 > m-^ , (18) 

will only be satisfied if the product of the quantum numbers N and Z is sufficiently large. A 
loop that does not satisfy this requirement will never stabilise as a vorton. After its length 
has been reduced to the order of magnitude ( [T8| ) by a classical contraction process, it will 
presumably undergo a rapid quantum decay whereby it will finally disappear completely 
just as if there were no current. 

III. THE VORTON ABUNDANCE 
A. Basic postulates: a scheme based on two mass scales. 

The present analysis will be carried out within the framework of the usual FRW model 
in which the universe evolves in approximate thermal equilibrium with a cosmological back- 
ground temperature T. The effective number of massless degrees of freedom at temperature 
T is denoted by g* . Note that (7* ^ 1 at low temperatures but that in the range where vor- 
ton production is likely to occur, from the electroweak scale through to grand unification, 
(7* ~ 10^ is a reasonable estimate. 

Any vorton formation processes must occur during the radiation dominated era which 
ended when the temperature of the universe dropped below 10^^ GeV and became effectively 
transparent. The relevant cosmological quantities are the age of the universe, given by 
if: H^^ where H is the Hubble parameter, and the radiation dominated time-temperature 
relationship 



where mp is the Planck mass. 

During this cosmological evolution, the particle physics gauge group is assumed to un- 
dergo a series of successive spontaneous symmetry breaking phase transitions, expressible 
schematically as 

g^,,^---n---^g,^^su{3)xu{i) . (20) 

Here ^qut "grand unified" group, Ti. is some hypothetical intermediate symmetry group 

(such as that of the axion phase), and Q^^^ is the standard model group SU{3)xSU{2)xU{l) 
or one of its non-standard (e.g. supersymmetric) extensions. The role of the Witten model 
is to provide an approximate description of an evolution process dominated by two distinct 
steps in this chain. 

When a semi-simple symmetry group, Q say, is broken down to a subgroup, 7i say, the 
topological criterion for cosmic string formation is that the first homotopy group of the 
quotient should be non trivial: 

MG/n} ^ 1 . (21) 

Our primary supposition is that such a process occurs at some particular cosmological tem- 
perature, Tx, which we assume to be of the same order of magnitude as the relevant Kibble 
mass scale m^. This mass scale is interpretable as being of the order of the mass of the Higgs 
particle responsible for the symmetry breaking according to the simple model discussed in 
the previous section. 

Our next basic postulate is that a current carrying field, characterised by the independent 
mass scale rria, condenses on the ensuing string defect at a subsequent stage, when the 
background temperature has dropped to a lower value, T^, which we assume to have the 
same order of magnitude as the mass scale m^. 

The formation of a condensate with finite amplitude characterised by the dimensionless 
sectional integral S does not in itself imply a non zero expectation value for the corresponding 
local current vector, j. However, one expects that thermal fiuctuations will give rise to a 
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non zero value for its squared magnitude |jp and hence that a random walk process will 
result in a spectrum of finite values for the corresponding string loop quantum numbers 

and Z . Therefore, in the long run, those loops for which these numbers satisfy the 
minimum length condition (0) are predestined to become stationary vortons, provided of 
course that the quantum numbers are strictly conserved during the subsequent motion, a 
requirement whose validity depends on the condition that string crossing processes later on 
are statistically negligible. We describe these loops as protovortons. 

Note that the protovortons will not become vortons in the strict sense until a lower 
temperature, the vorton relaxation temperature say (whose value will not be relevant for 
our present purpose) since the loops must first lose their excess energy. Whereas frictional 
drag and electromagnetic radiation losses will commonly ensure rapid relaxation, there may 
be cases in which the only losses are due to the much weaker mechanism of gravitational 
radiation. 

As the string network evolves, the distribution rarifies due to damping out of its fine 
structure first by friction and later by radiation reaction. However, not all of its lost energy 
goes directly into the corresponding frictional heating of the background or emitted radia- 
tion. There will always be a certain fraction, e say, that goes into loops which evolve without 
subsequent collisions with the main string distribution. It is this process that provides the 
raw material for vorton production. Such loops will ultimately be able to survive as vortons 
if the current induced by random fluctuations during the carrier condensation process is 
sufficient for the condition (^) to be satisfied, i.e. provided its winding number and particle 
number are large enough to satisfy 

\NZ\^I^ > ^ . (22) 

Any loop that fails to satisfy this condition is doomed to lose all its energy and disappear. 

In favorable circumstances, (namely those considered in Section IIIB) most of the loops 
that emerge in this way at will satisfy the condition (12^) and thus be describable as 



protovortons. However in other cases (namely those considered in Section IIIC) the majority 
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of the loops that emerge during the period immediately following the carrier condensation 
will be too small to have aquired sufficiently large quantum numbers by this stochastic 
mechanism. These loops will therefore not be viable in the long run and are classified 
as doomed loops. Nevertheless, even in such unfavourable circumstances, the monotonic 
increase of the damping length Lmin will ensure that at a lower temperature Tf < T„ a later, 
and less prolific, generation of emerging loops will after all be able to qualify as protovortons. 
We refer to Tf as the protovorton formation temperature. 

The scenario summarised above is based on the accepted understanding of the Kibble 
mechanism P] , according to which, after the temperature has dropped below the effect of 
various damping mechanisms will remove most of the structure below an effective smoothing 
length, Liniri) which will increase monotonically as a function of time, so that nearly all the 
surviving loops will be have a length L = § di that satisfies the inequality 

L ^ -^min (23) 

There will thus be a distribution of string loops, of which the most numerous will be relatively 
short ones, with L ^ -Lminj that are on the verge of emerging, or that have already emerged, 
as protovortons or doomed loops as the case may be. 

Whereas on larger scales closed loops and wiggles on very long string segments will be 
tangled together, on the shortest scales, characterised by the lower cutoff L^in, loops will 
be of a relatively smooth form. It is these smallest loops that are candidates for subsequent 
transformation into vortons. 

The total number density of small loops with length and radial extension of the order of 
Lmin will (due to the rapid fall off of the spectrum that is expected for larger scales) be not 
much less than the number density of all closed loops and so will be given by an expression 
of the form 

n^u L-f, (24) 
where i/ is a time-dependent parameter which we will discuss later. 
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The theory reviewed above was originally developed on the assumption that the string 
evolution is governed by Goto-Nambu type dynamics. In the kind of scenario we are con- 
sidering, this condition will obviously be satisfied as long as the cosmological temperature 
T is greater than or comparable with the carrier condensation temperature To-. Moreover, 
the usual Goto-Nambu type description, and its consequences as described above, will re- 
main valid for a while after the strings have become superconducting since the currents will 
initially be too weak to have significant dynamical effects. The Goto-Nambu theory in- 
evitably breaks down at some temperature above Tj.. However, there may be cases for which 
such a description will break down even before the protovorton formation temperature Tf is 
reached. 

The typical length scale of string loops at the transition temperature, L^in(T„), is con- 
siderably greater than relevant thermal correlation length, T~^, that will presumably char- 
acterise the local current fluctuations at that time. It is because of this that string loop 
evolution is modifled after current carrier condensation. The inequality 



will satisfy L > ^^111(^0-), means that the random walk effect can build up reasonably large, 
and typically comparable initial values of the quantum numbers \Z\ and A^. The reason is 
that for a loop of length L, the expected root mean square values produced in this way from 
carrier fleld fluctuations of wavelength A can be estimated as 



(25) 



and the fact that, by (p3D, the length of any loop present at the time of the condensation 




(26) 



At the time of the condensation, a typical loop is characterised by 




(27) 



and 



1 



(28) 
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so that one obtains the estimate 

\Z\^N^^L~(T^„, (29) 

which, by (^), is large compared with unity. 

For current condensation during the friction dominated regime discussed in Section IIIB, 
we shall see that this will always be sufficient to satisfy the requirement (^). However, this 
condition will not hold for condensation later on in the radiation damping regime discussed 
in Section IIIC. In the latter case, typical small loops that free themselves from the main 
string distribution at or soon after the time of current condensation will be doomed loops 
since they do not satisfy (^21) . However, there will always be a minority of longer loops for 
which (p^) will be satisfied, namely those exceeding a minimum length given, according to 
(H), by 

(30) 

(T 

This condition is still not quite sufficient to qualify them as protovortons since such excep- 
tionally long loops will be very wiggly and collision prone. It is not until a later time at 
a lower temperature Tf that free protovorton loops will emerge. In this case, the typical 
wavelength of the carrier field will be given by (p8|) and the final value of the length of a 
typical loop is 

L ^ L^UTi) . (31) 
The new estimate for the values of the quantum numbers is 



|Z|^iV^^ ^"^'°J'^^^'^ , (32) 

where we have included a blueshift factor, Z{, whose value is not immediately obvious 
but that is needed allow for the net effect on the string of weak stretching due to the 
cosmological expansion and stronger shrinking due to wiggle damping during the period as 
the temperature cools from T„ to Tf. In the earlier friction dominated regime Tf is identifiable 
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with Tcr so the problem does not arise, and in the radiation damping era cosmological 
stretching will in fact be negligible. Therefore, the net effect is that Zf will be small compared 
with unity, the hard part of the problem being to estimate how much so. 

The value given by (|32| ) will increase monotonically as Tf diminishes. The required value 
of Tf, at which the formation of the protovorton loops will actually occur, is that for which 
the function in (^) reaches the minimum qualifying value given by (^). This value is thus 
obtainable in principle by solving the equation 



(33) 



but this can only be done in practice when we have found the Tf-dependence of Ljnin(7f) 
and Zf. We discuss the Tf-dependence of Z^ shortly. 

The number density of protovorton loops at the temperature Tf will be comparable with 



the total loop number density at the time, so that by ( p4D it will be expressible as 



n,^ep,L^,{T,)-\ (34) 

where e is an efficiency factor of order unity, and z/f is the value of the dimensionless param- 
eter V at that time. If the current condenses in the friction dominated regime will simply 
have an order of unity value. However, if the condensation does not occur until later on, in 
the radiation dominated era, will have a lower value which is not so easy to evaluate. 

The number of protovorton loops in a comoving volume will be approximately conserved 
during their subsequent evolution. Therefore, since volumes will scale proportionally to the 
inverse of the entropy density it follows that the number density of the resulting vortons 
at a lower temperature T will be given in terms of the number density TZf of the proto-vorton 
loops at the time of condensation by 

=^-/f^r (35) 



Hi VTf. 

Here / is a dimensionless adjustment factor that we expect to be small but not very small 
compared with unity, and that will be given by 
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where is the value of g* at the protovorton formation temperature Tf. 



Using ( p!7D the corresponding mass density will be given by 



(37) 



Thus, the mass density of the distribution of the protovortons in the range Tf > T > Tf, 
and of the mature vortons after their formation in the range T < Tj., is given by the general 
formula 



When the dependences of z/f, Linin(Tf) and Tf on the fundamental parameters Tx and 
To- are known, the formula (|38| ) will allow us to place limits on by determining how 
the presence of the corresponding population of remnant vortons would affect the course 
of cosmic evolution. In the following sections we shall derive several constraints on such 
a population by demanding that it not significantly interfere with the cornerstones of the 
standard cosmology. However, before we can do so it remains to obtain at least rough 
estimates of the required values of the dependent variables. This turns out to be fairly 
easy in the case of condensation during the friction dominated era that will be discussed in 
the next section. However the derivation of firm conclusions is less straightforward for the 
kind of scenario discussed in Section IIIC, in which the current condensation occurs in the 
radiation dominated regime. 



tially dominated by the frictional drag of the thermal background. The relevant dynamical 
damping timescale, r, during this period is approximately given by 




(38) 



B. Condensation in the friction damping regime. 



According to the standard picture [^, the evolution of a cosmic string network is ini- 



(39) 
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where /3 is a dimensionless drag coefficient that depends on the details of the underlying 
field theory but that is typically expected p^ , [23[] to be of order unity. In this regime the 
large scale structure is frozen and retains the Brownian random walk form (p^. However, 
the microstructure is smoothed out below a correlation length Lmin given by 



^min~Vrt, (40) 

where t is the Hubble time. Neglecting the very weak (yf*-dependence, the required correlation 
length is thus found to be given by 

(41) 

The friction-dominated regime continues until the temperature T drops below a critical 
value given by 

^ TT^ , (42) 

at which r is comparable with t. 

Setting T equal to in ( |il| ) in order to obtain the relevant value of Lmin(To-), and using 
(PPP and (P7D, the required expectation value for the quantum number N can be estimated 

as 

/ \ /rp , 1/2 



It follows from (|17D and ([T6|) that a typical vorton in this relic distribution will have a 
mass-energy given by 



which corresponds to a vorton circumference 



^^'^^T.T.)-^/^ (45) 
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It can thus be confirmed using (^) that the postulate T„ > automatically ensures 
that these vortons will indeed satisfy the minimum length requirement ([I8|), though only 
marginally when T is at the lower end of this range. 

From ( p4D and (^) the number density of these proto-vorton loops, at formation, is 

It follows from (^) that at later times the number density of their mature vorton successors 
will be 

Thus, after the temperature has fallen below the value Tr, the resulting mass density of the 
relic vorton population will be 



p,^^JN\^\ (^)'t,T' , (48) 



-I VT 



which by gives our final estimate as 



C. Condensation in the radiation damping regime. 

For strings formed at low energies, for example in some non-standard electroweak sym- 
metry breaking transition, the scenario of the preceding subsection is the only one that needs 
to be considered. However, for strings formed at much higher energies, in particular for the 
commonly considered case of GUT strings, current condensation could occur during the 
extensive temperature range below Tj,. The minimum length requirement is only marginally 
satisfied by typical loops when condensation occurs near the end of the friction dominated 
regime. Therefore, if < T^, typical loops present during the transition will not be long 
enough to qualify as protovortons. This means that the vorton formation temperature Tf 
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will not coincide with T^- as it did in the friction dominated regime, but rather will have a 
distinctly lower value. 

In these scenarios the final stage of protovorton formation will be preceeded by a period 
of evolution in the temperature range > T > Tf. During this interval friction will be 
negligible and the only significant dissipation mechanism will be that of radiation reaction. 
Moreover, during the first part of this period, in the range > T > T^, the only radiation 
mechanism will be gravitational, which is so weak that to begin with it will have no percep- 
tible effect at all. Thus, there will be an interval during which the smoothing length remains 
roughly constant at its value at the end of the friction dominated era, given, according to 
(|1|) and (H), by 

Lr^UT.) ^ ^ . (50) 

X 

During the last stage before the protovortons are formed, in the range > T > Tf 
there will already be currents on the strings. However, in practice, even in the coupled 
case the expected currents will be too weak for electromagnetic radiation damping to be 
important. Therefore, gravitational radiation is the only important effect throughout the 
range > T > Tf. 

The resulting gravitational smoothing scale will be the length of the shortest loop for 
which the survival time exceeds the cosmological timescale (p^Qf). From dimensional consid- 
erations this can be estimated by the expression 

where £^ ~ T ~ is the mass energy density of the string. Here T is a dimensionless 
coefficient of order unity and, for the GUT strings the gravitational factor will be given by 
~ (mx/mp)^ ^ 10~^. The validity of the formula ( |5TD has been confirmed in many 
particular cases by numerical simulations 0, though the value of the coefficient turns out 
to be typically T ^ 10^. 

Equating (^) to the cosmological timescale we have 
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This formula is valid when its value becomes larger than that given by (0). This occurs at 



a critical value T| < Tj, which, from (52), is given by 




The relation (|52D may also be expressed as 



L 



mm 



(54) 



where k is a constant given by 




(55) 



To avoid ambiguity we can of course simply use the formula (^4|) as a defining relation 
to specify the parameter k during the Hubble damping "doldrum" regime > T > T'^. 
However, with this convention k will have to be considered as a function of time, starting 
with unit value at and decreasing to the very low value (|55D at which it levels off at 
T K,Ty This description of the network evolution is illustrated graphically in figure 1. 

In order to apply the formula (|38| ) for the final vorton mass density we must evaluate the 
dimensionless coefficient v that determines the protovorton number density. It was easy to 
do this for the friction dominated case for which we could assume a constant value of the 
order of unity. However, v may subsequently be reduced by an amount whose estimation is 
not so clearly evident. It is reasonable to expect that in the radiation dominated regime the 
string distribution tends towards a scale invariant form, albeit one that will not be quite so 
simple as the Brownian form that prevailed in the friction-dominated era. Although there 
may already be approximate scaling for large scales, it is clear that scaling on the smaller 
scales that matter for our present purpose can only be obtained after the parameter k defined 
by (|54D has settled down to a constant value. This occurs at temperatures below T-|- for which 
a detailed analysis is beyond the scope of the simulations that have been achieved so far. 
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For simplicity we shall use a crude, but we hope sufficiently robust, description based on 
simple and quite natural physical considerations. 

In so far as v is concerned, the plausible conjecture that it should be describable by 
a scaling solution is to be interpreted as meaning that it should be a function only of 
the dimensionless ratio Rjt (where, it is to be recalled, R denotes the radial scale under 
consideration and t is the Hubble time). It is reasonable to expect that for values of R in 
the range k <^ Rjt <^ 1, the value of v should be given by a simple power law of the form 



(56) 



t 

with constant index C,. We expect that, in the radiation dominated regime, the appropriate 
value of the index should be close to but perhaps slightly greater than a lower limit given 

by a 

C^f. (57) 

(The analogue for the matter era in which we are situated today is a value slightly greater 
than a lower limit given by C = 2.) 

Assuming that the formula ( [56| ) still gives the right order of magnitude at the lower end 
of its range, v will be given in the radiation dominated regime by the constant value 

z/ «i z/^ (58) 

with K given by (^Sp. This means that the corresponding value of the loop number density 



itself will be given according to ([5^ ) by 

nK,v^ K^-H'^ . (59) 

Before evaluating the required result (pS]), it remains to obtain the value Tf. To do this 
we have to solve the equation (^Sj) that results from the minimum length requirement, which, 
from (|52|), reduces to 

Z,T,' ^ — =^ . (60) 
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However, before we can solve this deceptively simple equation in practice, we need to know 
the T dependence of the factor Z. This is the most delicate part of the calculation, since it 
involves competing effects of comparable magnitude. 

As the string distribution evolves, the time dependent blue shift factor Z is the factor 
by which the supporting string length has shrunk since the time of the condensation at the 
temperature To-. This shrinking can be accounted for as the net result of three main effects, 
two of which are comparitively easy to evaluate. 

The weakest of these effects is the stretching due to the expansion of the universe. 
This will always be more than compensated (except in the "doldrum" period in which 
compensation is only marginal) by shrinking due to the steady damping out of the short 
wavelength wiggles that give the most important contribution to the total string length per 
unit volume in the radiation dominated era. If these two effects were the only ones it would 
be relatively easy to estimate Z since it would simply be proportional to the total string 
length, H say, in a comoving volume that can conveniently be taken to be a cubic thermal 
wavelength. This is given by 

E ^ AT-3 (61) 

where A is the total string length per unit volume. Now, note that the main contribution to 
the total length of the string distribution is provided by short wavelength modes with scale 
of order the smoothing length L^in- Therefore, A can be estimated as being nLmin, where n 
is given by (p^ . This implies 

A^z/L-2^, (62) 

in which u is given by (0). 

If the only effect of the damping were to smooth out the short wavelength wiggles on 
the main part of the string distribution, it would be possible to identify Z with the ratio of 
S to its value S^j at the time of the current condensation. However, this would not allow 
for a third important effect, namely the losses of small loops which are continually liberated 
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from the string distribution at the lower end of the spectrum. Whether these loops survive 
as vortons or disappear altogether, the effect on the main part of the string distribution 
is that the corresponding string lengths must be subtracted at each stage. Thus, the total 
shrinking factor Z will be the ratio of Z", not to its original known value but to a value 
that is considerably reduced in such a way as to take account of this new effect. 
In terms of the variation of 17, the variation bZ of Z is given by 

bz (5i: + Ar 

where AZ" is the length of string irreversibly chopped off in the form of small loops per 
comoving thermal volume within the short time interval bt under consideration. The delicate 
question is that of quantifying AZ. It is not hard to obtain an order of magnitude estimate, 
but since the final result is rather sensitively dependent on this quantity a more accurate 
estimate would be desirable. In the absence of a detailed numerical investigation, we express 
AZ as a fraction of the total lost length 

AZ ~ -ebE , (64) 

where e is a dimensionless efficiency factor that must lie in the range < £ < 1. This 
factor is roughly identifiable with the coefficient introduced, using the same notation, in 
(0). However in that context it was sufficient to know that it should be of order unity. 

Whatever the exact value of e, the substitution of (0) in (|63|) provides a differential 
equation that can be solved to give 

Z (65) 



and using (p2|) and (p^) we finally obtain 



^-(^)"^ (66) 

It is to be remarked that if the efficiency e of loop production were zero, this would mean that 
the carrier field would be blue shifted by a factor that would be precisely the inverse of that 
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by which the background radiation is redshifted. However, if substantial loop production 
occurs there will be a blueshift by a moderate factor. 

Assuming the particle number weighting factor can be taken to have the fixed value 
g*, ( |60D and (p6| ) then give 



Using this result in conjunction with the estimates and the formula gives 
the mass density of the resulting vorton distribution as 

p. « ,^v.r^-/^(^n^-'Vl^)='=(^)"«'"=-"' rj' . (68) 



If we adopt the value (|57D for this simplifies to 

but there remains an unsatisfactory degree of sensitivity to the uncertain efficiency factor 
e. This is because although this index will always be quite small, the factor T„/mp is very 
tiny in the cases of interest. 

D. Stability Issues 

Before we consider cosmological constraints we would like to say a few words about 
stability. One of our postulates is that superconducting current conservation is sufficiently 
effective to allow the protovortons that emerge at the temperature Tf < T^- to settle down 
as dynamically stable bodies, at a possibly lower relaxation temperature ^ Tf. This does 
not exclude the possibility that in the very long run they may finally decay by quantum 
tunneling or other "secular" instability mechanisms. In this case they would ultimately 
disappear when the thermal background reached a corresponding vorton death temperature 
with an even lower value, Td say. This means that a complete analysis of vorton formation 
and evolution could involve five successive temperature scales related by 

Tx > > Tf > T, > Td . (70) 
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Protovortons are small compared with the ever expanding scales characterising the rest of 
the string distribution and hence undergo few extrinsic collisions. Also, they are sufficiently 
smooth to avoid destructive fragmentation by self collisions. It therefore follows that in 
most cases the relevant quantum numbers N and Z will be conserved. As a consequence, 
the statistical properties of the future vorton population will be predetermined by those 
of the corresponding protovorton loops at the time of their emergence at the temperature 
7f. It is therefore unnecessary for our present purpose to consider how long it takes for the 
protovorton loops to settle down and become proper stationary vortons. Thus, the value 
of Tr will not play any role in the discussion that follows. This is convenient because the 
details of protovorton loop decay have not yet been adequately studied, and will obviously 
be sensitively dependent on whether the current is electromagnetically coupled, in which 
case would expect the later stages of the protovorton loop contraction to be relatively rapid. 

The final decay temperature is also absent from the quantitative formulae that we 
will derive. Indeed, its only role is to characterise the two principle kinds of scenario that 
we consider. In Section V, we discuss vortons which survive until the present epoch, which 
requires that should not much exceed 10~^^ GeV (corresponding to the observed 3 degree 
radiation temperature). However, in Section IV we adopt the weaker supposition that the 
vortons survive at least till the time of nucleosynthesis, which requires only that should 
not much exceed 10~^ GeV. This means that the only temperature scales that remain as 
variable parameters in the analysis that follows are the Higgs-Kibble temperature Tx, the 
condensation temperature and the protovorton formation temperature Tf which will not 
be truly independent but is a function of 7^ and Tx. 

IV. THE NUCLEOSYNTHESIS CONSTRAINT. 

One of the most robust predictions of the standard cosmological model is the abundances 
of the light elements that were fabricated during primordial nucleosynthesis at a temperature 
^ 10-^ GeV. 
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In order to preserve this well established picture, it is necessary that the energy density in 
vortons at that time, pv{T^) should have been small compared with the background energy 
density in radiation, ^ 9*T^- Assuming that carrier condensation occurs during the 
friction damping regime and that g* has dropped to a value of order unity by the time of 
nucleosynthesis, it can be seen from (HUl) that this restriction. 



Pv(Tj « , (71) 

is expressible as 

eu.g:-'p'/'mp-'/'T-'/X''/' « . (72) 
Below we apply this constraint to some specific examples. 

A. Case Tx « T^. 

The case for which carrier condensation occurs at or very soon after the time of string 
formation has been studied previously and yields rather strong restrictions for very long 
lived vortons 0. If it is only assumed that the vortons survive for a few minutes, which is 
all that is needed to reach the nucleosynthesis epoch we obtain a much weaker restriction. 
Setting To- equal to in (^2]) gives 

/ \ 4/9 / s 5/9 

(t) --<(?) 

Taking g* ^ 10^ and assuming (in view of the low value of the index) that the net efficiency 
factor (ez/^)^/^ and the drag factor are of the order of unity yields the inequality 

Tx < 10^ GeV . (74) 

This is the condition that must be satisfied by the formation temperature of cosmic strings 
that become superconducting immediately, subject to the rather conservative assumption 
that the resulting vortons last for at least a few minutes. It is to be observed that this 
condition rules out the formation of such strings during any conceivable GUT transition, 
but is consistent by a wide margin with their formation at temperatures close to that of the 
electroweak symmetry breaking transition. 

26 



B. Case Tx 



=i T, 



GUT 



10^^ GeV. 



Here we wish to calculate the highest temperature at which GUT strings can become 
superconducting without violating the nucleosynthesis constraints. Setting equal to T^^^ 
in ([72|) , and again using g* ^ 10^, we obtain 



where, in the last step, we have neglected the dependence on order of unity quantities. It 
can be checked, using the Kibble formula (^2]), that the maximum value given by ([75|) is 
at least marginally consistent with the assumption that current condensation occurs in the 
friction-dominated regime. The validity of our derivation is thereby confirmed. It follows 
that the nucleosynthesis constraint will always be satisfied when T^. lies in the radiation 
damping epoch. 

Therefore, subject again to the rather conservative assumption that the resulting vortons 
last for at least a few minutes, theories in which GUT cosmic strings become superconducting 
above 10^^ GeV are inconsistent with the observational data. 



In this section we consider the rather stronger constraints that can be obtained if at 
least a substantial fraction of the vortons are sufficiently stable to last until the present 
epoch. It is generally accepted that the virial equilibrium of galaxies and particularly of 
clusters of galaxies requires the existence of a cosmological distribution of "dark" matter. 
This matter must have a density considerably in excess of the baryonic matter density, 
Pb ~ 10~^^ gm/cm'^. On the other hand, on the same basis, it is also generally accepted that 
to be consistent with the formation of structures such as galaxies it is necessary that the 
total amount of this "dark" matter should not greatly exceed the critical closure density, 
namely 




(75) 



V. THE DARK MATTER CONSTRAINT. 



Pc ~ 10 gm cm 



-3 



(76) 
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As a function of temperature, the critical density scales like the entropy density so that it 
is given by 

p,(T) ^ g*m,T^ , (77) 

where mc is a constant mass factor. Since g* 1 aX the present epoch, the required value 
of TTic (which is roughly as the critical mass per black body photon) can be estimated as 

mc ^ lO^^^mp « 1 eV . (78) 

However, for comparison with the density of vortons that were formed as a result of current 
condensation at an earlier epoch characterised by T^., what one needs is the corresponding 
factor g*mc, which can be estimated to be 

glrric ^ lO^^^mp ^ 10^ eV . (79) 

(This distinction was obscured in the previous derivations of the dark matter constraint 
pO| , in which the value quoted for nic should be interpreted as meaning the value of g*mc, 
which is what one actually needs.) 

The general dark matter constraint is 

f2, = ^<l. (80) 

Pc 

In the case of vortons formed as a result of condensation during the friction damping 
regime the relevant estimate for the vortonic dark matter fraction is obtainable from (^Of) as 

In particular, this formula applies to the case in which the carrier condensation occurs very 
soon after the strings themselves are formed, as was supposed in earlier work. 

However, if we want to strengthen the nucleosynthesis limit (^) for the general category 
of strings formed at the GUT scale, then we are obliged to consider the case of vortons 
formed as a result of condensation during the gravitational radiation damping regimes. In 
this case, equation (p9| ) gives the relevant estimate for the vortonic dark matter fraction as 

Let us now again examine some specific examples. 
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A. Case T^^T^. 



The formula (^Tj) is applicable to the case considered in earlier work , in which it was 
supposed that vortons sufficiently stable to last until the present epoch were formed as the 
result of the carrier condensation occurring at or very soon after the time of string formation. 
This example provides the strongest limits on T^. Setting T„ equal to in (|ST| ) one obtains 




Substituting the estimates above (supposing, as before, that the efficiency and drag factors 
are order unity), we obtain 

< 10^ GeV . (84) 

This result is based on the assumptions that the vortons in question are stable enough 
to survive until the present day. Thus, this constraint is naturally more severe than its 
analogue in the previous section. It is to be remarked that vortons produced in a phase 
transition occurring at or near the limit that has just been derived would give a significant 
contribution to the elusive dark matter in the universe. However, if they were produced at 
the electroweak scale, i.e. with Tx ~ T^. ~ where T^^ ^ 10^ GeV, then they would 
constitute such a small dark matter fraction, f]^ ^ 10"^, that they would be very difficult 
to detect. 

B. Case ~ T^^^ ^ 10^^ GeV. 

For the most commonly considered case, namely that of strings formed during the GUT 
transition, the nucleosynthesis limit ([75| ) is already sufficient for the exclusion of carrier 
condensation in the friction damping regime. To obtain the stronger limit that is applicable 
if the vortons are sufficiently stable to survive as a dark matter constituent, we need to 
consider the case in which the condensation occurs during the the regime of gravitational 
damping. In this case, the relevant dark matter fraction is given by (^). Setting Tx equal to 



29 



T^yrj, in this formula, and dropping the order of unity coefficients we obtain the corresponding 
hmit 



^ < /mcT^N(3-£)/(9-2e) 

If the loop production were extremely efficient, e ~ 1, this would already give the numerical 
limit 

< 10^° GeV . (86) 

which is significantly stronger than the more conservative limit ([75| ) that pertains if the 
vortons only survive for a few minutes. 

However, contrary to what one might have guessed, the highest conceivable loop produc- 
tion efficiency is not what maximises ultimate vorton production. This is because it merely 
tends to enhance the charge and current loss rate by production not of protovortons but of 
doomed loops. Thus if, instead of supposing that the loop production efficiency e is close to 
a hundred percent, one makes the plausible supposition that it does not much exceed fifty 
per cent, then one obtains 

< 10^ GeV . (87) 

This limit is still compatible by a very large margin with one of the most obvious, albeit 
rather extreme, possibilities ^ that comes to mind, namely that in which the strings are 
formed at the GUT level, Tx ~ T'gut; but no current carrier condenses until the electroweak 
level, ^ T^^. 

VI. CONCLUSIONS 

We have explored the constraints and implications both for particle physics and cos- 
mology arising from the existence of populations of remnant vortons of more general types 
than have previously been considered. Specifically, we have envisaged the possibility of cos- 
mic string superconductivity by condensation of the relevant carrier field at energy scales 
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significantly below that of the string formation. We have seen that there are two quali- 
tatively very different possibilities. In scenarios for which the carrier condensation occurs 
at comparitively high energy, during the friction damping regime, a substantial majority 
of the superconducting string loops will ultimately survive as vortons. Such scenarios are 
more easily excluded on observational grounds than the alternative possibility, which is that 
superconductivity does not set in until a later stage, in which case only a minority of the 
loops initially present ultimately become vortons. 

We have shown that large classes of particle physics models can provisionally be ruled out 
as incompatible by these cosmological considerations, and in particular we have shown that 
models admitting GUT strings must not allow any string superconductivity giving stable 
vortons to set in much above 10^ GeV. The excluded regions of parameter space are shown 
in figure 2. 

Our conclusions are, however, dependent on a number of more or less "conventional" 
assumptions, whose validity will need to be systematically scrutinised in future work. In- 
validation of these conventional assumptions, particularly those concerning the long term 
stability of the vortons, in specific theoretical contexts would mean that in such circum- 
stances the constraints given here might need to be considerably relaxed. On the other 
hand, our constraints may be considerably tightened by the use of more detailed observa- 
tional data and the ensuing limits on the populations of various kinds of vortons that can 
exist today. On the constructive side, we have shown that it is possible for various conceiv- 
able symmetry breaking schemes to give rise to a remnant vorton density sufficient to make 
up a significant portion of the dark matter in the universe. 
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Figure Captions 

Figure 1: The evolution of relevant physical quantities as a function of the thermal length 
scale in the early universe. 

Figure 2: A summary of our results. This diagram shows how the appropriate length scales in 
particle physics theories are constrained. Both the nucleosynthesis and dark matter bounds 
are shown. 
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Vorton formation as a function of superconducting 
transition temperature for GUT strings. 
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